Josephson current between chiral superconductors 

Yu. S. Barash 1 - 2 and A. M. Bobkov 2 

1 Center for Electronic Correlations and Magnetism, Institute of Physics, University of Augsburg, D-86135 Augsburg, Germany 

2 P.N. Lebedev Physical Institute, Leninsky Prospect 53, Moscow 119991, Russia 



o 

o 

o 
o 

or 

+-> 
c3 



-a 

c 
o 

> 

in 

o 
l> 

o 



C3 



I 

c 
o 
o 



s 



M. Fogelstrom 
Institute of Theoretical Physics, Chalmers University of Technology and Goteborg University, S-J i 1296 Goteborg, Sweden 

We study chiral interface Andreev bound states and their influence on the Josephson current 
between clean superconductors. Possible examples are superconducting Sr2Ru04 and the B-phase 
of the heavy-fermion superconductor UPt3. We show that, under certain conditions, the low-energy 
chiral surface states enhance the critical current of symmetric tunnel junctions at low temperatures. 
The enhancement is substantially more pronounced in quantum point contacts. In classical junctions 
dispersive chiral states result in a logarithmic dependence of the critical current. This logarithmic 
behavior contains the temperature, the barrier transparency and the broadening of the bound states, 
and depends on the detailed relation between these parameters. The Josephson current through the 
domain wall doesn't acquire this logarithmic enhancement, although the contribution from the 
bound states is important in this case as well. 

PACS numbers: 74.50. +r, 74.80. Fp 



I. INTRODUCTION 

Superconductors in a state with both time-reversal 
symmetry and parity with respect to the inversion of 
a crystal axis broken are sometimes named chiral su- 
perconductors. The "chirality" can be associated with, 
for example, a nonzero projection of the orbital angu- 
lar momentum, I, of Cooper pairs along a axis z, i. e. 
l z 7^ 0. Superconducting Sr 2 Ru04 as well as the heavy- 
fcrmion superconductor UPt3 in its i?-phase (i.e. the 
low-temperature and low-field phase) are leading candi- 
dates for chiral superconductors of this type. Another 
possible reason for a nonzero chirality of a superconduct- 
ing state to appear is a complex admixture of two pair- 
ing channels belonging to different irreducible representa- 
tions. For two-dimensional chiral superconductors there 
is a topological invariant associated with the chiralityEl. 

One of the important features of superconductors is the 
presence of surface or interface Andreev bound states. 
They can take place in the vicinity of surfaces or in- 
terfaces if the order parameter varies in space or takes 
different values for incoming and outgoing momentum 
directions along a quasiclassical trajectory travelled by 
a quasiparticle in a reflection or a transmission event. 
Surface and interface quasiparticle states arising in su- 
perconductors already in the absence of the magnetic 
field, have been -studied theoretically for many years 
(for example, seeafl). They have attracted much atten- 
tion in investigating high temperature superconductivity. 
For instance, the zero-energy surface states in d-wave 
superconductors appear near an impenetrable smooth 
surface as a consequence of a sign change of the or- 
der parameter along quasiparticle trajectories connected 
by the scattering off the surface. Zero-energy bound 
states lead at low temperatures to a zero-bias conduc- 
tance peak of N-FD (normal metal - isolating barrier - 



el junctionsBLl, which was 
Low-temperature anoma- 



d-wave superconductor) ± 
observed in experimentsB 
lies, originating from the zero-energy bound-states, are 
also seen in the JosepJj^on-critical currentcj cj and in 
the penetration depthEj£3MEZI. The strength of these 
effects depend on the .broadening of the bound [States 
by surface roughnessEilEirLJ and bulk impuritiesEj, as 
well as on their shift from zero on account of junction 
transparencyEj'E3E3. 

Zero-energy surface states are dispersionless states, 
while Andreev bound states with nonzero energy are dis- 
persive: the energy depends on a quasiparticles momen- 
tum direction. Bound states at energies on the order 
of A(pf) can manifest themselves, for example, in the 
current- voltage characteristics, pi a junction along with 
the low-energy bound statesLrEl. They do not, how- 
ever, modify noticeably equilibrium characteristics like 
the Josephson critical current or the penetration depth. 
In contrast, the zero- and/or low-energy states can have 
strong influence on the equilibrium quantities at low tem- 
peratures, ni— ii— i 

Chiral surface statesErEilc3 make up a special type of 
Andreev bound states. They form dispersive branches, 
which can cross the zero level only for isolated quasi- 
particle momentum directions. The-sppntaneous sur- 
face current in chiral superconductorsal£3, the spectrum 
of thc-surface chiral bound states at an impenetrable 
wallu't 2 ], and their contributions to the density-of-states, 
the conductance peak of the tunnel jurptionsElTEil, and 
the low-temperature penetration deptbo, have already 
been studied theoretically in literature. There are, how- 
ever, various open problems in the field, part of which is 
addressed in the present paper. We report on theoreti- 
cal results for the spectra of chiral surface and interface 
states and their contribution to the Josephson critical 
current of junctions between chiral superconductors. In 



experiments to date the Josephson effect has been stud- 
ied only in junctions betweep-&eRu04 (or UPts) and a 
conventional superconductoreS'EZl. Our results show that 
extending experimental investigations of the Josephson 
critical current in symmetric junctions of the ruthenates 
(or UPta) would sensitively probe the chiral nature of 
these superconductors. 

Our numerical approach is to solve the quasiclassi- 
cal transport equation and to determine the profile of 
the order parameter self-consistently. Our analytical re- 
sults describe spectra of the chiral states at interfaces 
of any transparency assuming spatially constant (non- 
self consistent) order parameters. These calculations are 
improved further for low-energy states in the tunnel- 
junction limit, where the surface pair-breaking may be 
taken into account analytically as well. As is demon- 
strated below, the low-energy parts of the chiral branches 
of the surface Andreev bound state spectra can dom- 
inate and strengthen the Josephson critical current in 
symmetric tunnel junctions at sufficiently low tempera- 
tures. The low-temperature enhancement of the Joseph- 
son critical current, j c , is substantially more pronounced 
for quantum point contacts. It is of similar origin as the 
low-temperature anomaly, in j c of junctions between d- 
wave superconductorsE3Eil. In classical junctions the dis- 
persion of the chiral states strongly modifies the anoma- 
lous low-temperature behavior of j c , leading eventually 
to an additional logarithmic factor, which can be large 
under certain conditions: if the broadening of the bound 
states and the transparency of the junction are suffi- 
ciently small. The zero-temperature value of j c is then 
determined by a presumably small broadening of the 
bound states or/and a small junction transparency that 
cuts off the logarithmic divergence as T — > 0. 

We also study the Josephson current through a domain 
wall, considering a junction between identical equally 
oriented superconductors with opposite chiralities. We 
demonstrate the crucial role of the Andreev interface 
states for obtaining a finite Josephson current in the 
system with the crystal c-axis parallel to the interface. 
Disregarding the presence of surface states and assum- 
ing spatially constant order parameters, one can apply at 
the interface the bulk expression for the Green's function. 
Then the tunnel Josephson critical current across the do- 
main wall vanishes. If the c-axis is along the interface 
normal and the interface itself is symmetric with respect 
to rotations around the normal, then the dc Josephson 
current actually vanishes, which is known for a long time 
for junctions between superconductors with opposite p.ra- 
jections of the angular momentum of Cooper pairsl3E3. 
For the c-axis parallel to the surface, vanishing current is, 
however, a shortcoming of the oversimplified approach. 
We show, that the interface influence results in this case 
both in Andreev bound states and in a finite Josephson 
current close to the conventional value. 

We perform our analytical and numerical calculations, 
first, for a quasi two-dimensional tetragonal superconduc- 
tor, choosing several particular basis functions, which are 



considered as candidates for the superconducting pairing 
in S^RuOzi. Possible candidates for the order .pa ra meter 
in Sr2Ru04 are still intensively discussecEJE§I^E2~ii3. 
For our study we take three types of two-component 
triplet order parameters forming the p-wave, the f xy - 
wave and the f x 2_ y 2 superconducting (l,±i) phases. 

We also present analytical results, demonstrating 
that analogous chiral interface bound states and low- 
temperature enhancement in the Josephson current take 
place in a three-dimensional hexagonal superconductor 
like UPt3. For many years the two-component triplet 
^^-representation and singlet i?i 9 -representation order 
parameters were considered .as ±he probable candidates 
for the pairing state in UPt3E3~E3. In particular, (1, ±i) 
superconducting phases are supposed to form at low tem- 
peratures and in weak magnetic fields. The analysis of 
experimental data now available has led to a presum- 
ably definite, choice in favor of the £211 type of pairing 
in UPt 3 ErO. The so-called "3D" model for the order 
parameter in tetragonal Sr2Ru04 is based on an analo- 
gous type of pairingpl We examine pairing both in the 
E-iu and in the E\ g representation for three-dimensional 
superconductors and find conditions, when only one of 
them results in the anomaly of the Josephson current we 
study below. 



II. CHIRAL SURFACE AND INTERFACE 
QUASIPARTICLE STATES 

A. Chiral order parameters 

Let a superconductor breaking time-reversal symme- 
try occupy a half-space with a surface (or interface) nor- 
mal along the :r-axis. The superconductor may be both 
triplet or singlet. Triplet order parameters, A{pf) = 
d(pf)-& ia y , considered below have only the component 
d z {pf) = A(pf) nonzero. 

Unconventional order parameters are in general quite 
sensitive to any inhomogeneity in the superconduc- 
tor that gives rise to quasiparticle scattering and this 
in particular to a boundary. We represent a chi- 
ral order parameter near a surface or interface as 

A(p f ,x) = (Ai(p / ,:r)e^ 1 ^-MxA2(p/,x)e^ 2 ( x ))e l * = 

\A(pf,x)\e l ^ Pf ' x \ where the two real components 
A\(pf,x) and A2(p/,x) have presumably opposite par- 
ities with respect to the inversion of the x-axis. The 
parameter x = ±1 corresponds to two possible values of 
"chirality". The phases <f>i(x) and 4>2(x) vanish in the 
bulk but may be nonzero in vicinity of the boundary. 

Assuming a model with a single cylindrical Fermi sur- 
face, the two components of the p-wave order parame- 
ter in a tetragonal superconductor with the crystal xq- 
axis along the surface normal can be taken in the bulk 
as A p lb {pf) — A^cosy and A^ b {pf) — A^simp. Here 
and below tp is the angle which the vector {pf,x,Pf,yify 



makes with the x-axis. For the f xy -w&ve order parame- 

f f 

ter in the bulk, one has A^jf (pj) — A b xy sin 2<p cos (p and 



Jxy 



A0(P/) = A 
order parameter: 



sin 2ip sin ip. For the bulk f x 2 -y 



a(; & 2 - 2 (p/) 



cos 2ip cos ip 



and A 2 X 6 H (p/J = A^ x " y " cos2( / 9sin< / o 

If </3 describes the incoming momentum direction spec- 
ular quasiparticle reflection from the surface gives (tt— <p) 
for the outgoing momentum. This corresponds to the 
inversion of the x-axis. One can easily see that order 
parameter components A\(pf) and ^{Pf) have oppo- 
site parities with respect to the inversion of the x-axis 
in all cases considered. Since the total order parame- 
ter represents a complex mixture of A±(pf) and A2(p/), 
we get, indeed, chiral order parameters describing super- 
conducting states with both time-reversal symmetry and 
parity with respect to the inversion of the x-axis bro- 
ken. An additional feature, which makes our analytical 
results comparatively simple, is that in all cases we con- 



sider |A&(p/)| doesn't change in the reflection event (if 
01.2(20 = 0, then \A(pf,x)\ doesn't change). 

In Figure [fl we show order-parameter profiles at an 
impenetrable barrier separating a Ai(pf,x)e l< ^ 1 ^ + 
iA 2 {p f ,x)e lMx *> (left side) and a A 1 {p f ,x)e l(t ' 1 ^ - 
iA 2 {pf,x)e' l 't >2 ( x ' > (right side) superconductor. The tem- 
perature is 0.001 T c . The order parameter compo- 
nents are represented in the forms Ai(p/, x)e l ^ >1 ^ = 
Ai^Gv), ±iA 2 (p f ,x)e i ^ = A 2 M (z)V> 2 (P/), 
which define the profiles A : r 2 '{x) displayed in Fig. 
p]. They are computed using a single cylindrical 
Fermi surface and the pairing potential V(p, p') = 
Vo[ipi(p)ipi(p') + "02 (p) ^2(2?')] in a factorized form with 
the simple, normalized, basis functions 



(Mpf),Mp/)) = 



v2 (cos ip, sin ip) (p-wave) 

2 (cos</?, sinc/j) sin 2<p (f xy -wave) 



(1) 



ImA 



2 / 



0.2 



Re A, 



0.1 - 



0.0 



p-wave 



-0.1 



Im A 7 



Re A, 




Re A, 



ImA, 



f -wave 

•> xy 



Re A, 



ImA, 



Re A i ==Im A 2 




f 2 2 -wavfc 

J x -y 



-30 -20 -10 10 20 30-30 -20 -10 10 20 30-30 -20 -10 10 20 30 

d%o X ^o ^o 

FIG. 1. The spatial profiles of the order-parameter components Ai^z) (in units 2irT c ). 



For the p-wave case our self-consistent numerical re- 
sults are in agreement with those represented ir£3. One 
should note the vanishing Ai (x, pj) at the surface and, as 
a consequence, the increase in A2(0,p/) as compared to 
its bulk value. An inhomogeneous profile Ai(x,p/) near 
the surface slightly modifies the surface value A2(0,p/) 
as wellE3. If the crystal axis xo makes a finite angle a in 
the xj/-plane with the surface normal, for the p-wave su- 
perconductor there is no dependence on a. The surface 
state of an f-wave superconductor, however, will depend 
on a. The surface profiles of the /-wave order parameters 
with angles a — and a — 45° are shown in Fig. [H and 



correspond to the most extreme cases, f xy and jx^-y 2 re- 
spectively. We note, that the / xy -wave order parameter 
behaves quite similar to the p-wave case near a surface. 
For p-wave and f xy -w&ve pairings the phases <pi (x) and 
4>2(x) obey 4>i{x) = 4>2{x) = 0, while for the f x 2_ y 2-w&ve 
order parameter 4>\{x) = ~4>2{x) 7^ 0. For all pairing 
states we see that the order parameter has healed after 
roughly 15 £0, (£0 = v f/^ n T c ). If instead the right super- 
conductor would have been calculated with the same chi- 
rality as the left one, the order parameter would change 
as (A 1 (p f ,x),A 2 (pf,x)) -> (Al(p f ,x),-A*(p f ,x)) on 
the right side. When we continue to study the Josephson 



coupling of two chiral superconductors we shall assume 
that the Josephson currents generated are much smaller 
than the bulk critical current and do not affect the spa- 
tial dependence of the order parameter or the surface 
currents (see below). 

For three-dimensional models of chiral supercon- 
ducting states we assume, for simplicity, a spheri- 



cal Fermi surface and choose bulk order parameters 



in the form: 



A*, 



(pi) 



= A* cos 9 sin 2 Q e ^xf^) f or 



the ^u-representation (triplet pairing) and Ag(p/) = 
A£cos#sin#e^ xv+ *) for the E\ g (singlet pairing). 
Hence, for the Ei u and the E\ g pairings respectively we 
get 



A* b (pj) = A^ cos sin 2 6cos2tp , 
A* 2 b {pf) = A\ cos 9 sin 2 9s'm2p , 



Af b {pf) — Ag cos 9 sin 9 cos p , 
A 2 b {pf) — Al cos 9 sin 9 sin^p . 



(2) 



Here 9 and p are standard angles in a spherical system of coordinates in momentum space with the x-axis along 
the surface normal. $ is a spatially constant part of the complex phase of the order parameter. Order parameter 
components A|(p/,x) and Af(pf,x) vanish on an impenetrable surface. 



B. Interface chiral bound states 



Consider a junction between identical equally oriented chiral superconductors. Solving Eqs.(A4)-(A6) of the Ap- 
pendix with spatially constant anisotropic order parameters on both sides of the interface we find comparatively 
simple analytical expressions for the energies of chiral interface Andreev bound states. If x = X l X r = 1> i- e - the 
superconductors have identical chiralities, we get the following spectrum for two-dimensional p- wave and f x 2_ y 2-wave 
order parameters as well as for the three-dimensional singlet i?i g -pairing: 



_p,4 



(P/)=± 



A 



P.4 



For the E% u - and the /^-pairing we obtain 



(Pf] 



1Z sin 2 p + V cos 2 ( — 



e t B (p f ) = ±\At(p f )\ 



Tl cos^ 2p + V cos^ ( - 



1/2 



1/2 



(3) 



(4) 



>/) 



± 



'(Pf) 



11 cos^ p + V cos^ ( - 



1/2 



(5) 



Here T> and 1Z are transmission and reflection coefficients of the interface barrier which in general depend on the 
quasiparticle momentum direction. The phase difference <& = $ r — $', where ^ r ' 1 are the phases of the order 
parameters in the right and the left superconductors, is separated from the intrinsic chiral phases, ip, as defined in 
the preceding subsection. The spectra in Eqs.(|3|)-(|5|) depend on momentum directions both explicitly and via 23-and, 
1Z. They significantly differ from the bound state energies in junctions between conventional superconductorsE3"LJ 

e B {p f ) = ±\A\^Jl-Vsm 2 ^. 

Positive and negative branches in Eqs.(g)-(Q) always touch on the Fermi surface (on the zero-energy level) for 
momentum directions along possible nodes of |A(p/)|. In addition, if the phase difference <& is equal to 7r (or for any 
<f>, but in the limit of an impenetrable wall V = 0), the branches described by Eqs.([|) and (Q), acquire new crossing 
points (or lines) lying on the Fermi surface. These additional momentum directions are at <p = 0, it for the p-wave, 
the f x 2_ y 2-wave and the three-dimensional -Bi g -pairing, while at p — ±7r/4, ±37r/4 for the -E 2u -pairing. In accordance 
with Eq. (ffl) for the /^-pairing the spectrum of the bound states takes zero value only for directions along the nodes 
of the order parameter, i.e at ip = 0, 7r, ±7r/2. At the glancing trajectories p ~ ±7r/2 the nodes in the spectrum have 
a higher multiplicity. 

It is worth to compare interface states in Eqs.(|3|)-(|5|) with the respective surface chiral branches at an impenetrable 
wall: 






~ y2 (Pf) =xsgn(:r) 



p^J^-i. 



"(Pf) 



sin p 



(G) 



'(Pf) = -XSg'a{x)A b n:y sin 2ip cos 



Y 



Api) 



-XSgn{xp y )\A t b (pf)\cos2p> 



(7) 



One can see, that finite transmission and $ ^ 7r re- 
sult in an effective "repulsion" and a reconnection of 
the branches from two sides of the barrier plane. Chiral 
states on an impenetrable surface of the p-wave super- 
conductor were obtained earlier in Refs. plLp4 In the 
particular case $ = tt Eqs.(|g) and (J7|) describe as well 
the states on interfaces with arbitrary transparency if 
compensated by a factor yK. 

The chiral spectra in Eqs.(^|) and ([?]) are antisymmet- 
ric with respect to the inversion p y — > — p y (or equiva- 
lently ip — ► — (p), as well as to the inversion of the total 
momentum pf — > — p/ (or ip — > p+n). Since at zero tem- 
perature the states with negative (positive) energies are 
occupied (empty) the chital states carry a spontaneous 
current along the surfaceErE 2 ]. This can be considered as a 
characteristic feature of chiral states. The current flows 
near the surface parallel to the y-axis within a thickness 
the order of the coherence length £q. With the model of 
a constant order parameters we find at zero temperature: 
j y = — £%sgna; evjNfAb, where £ p — 1, (J^-v* = 0.5, 
£/** = -1, C, 1 « 0.08, C = 0.25. The induced screen- 
ing supercurrents decay on the order of the penetration 
depth, A, so that the total surface current is actually zero. 

As it follows from Eqs.(0) and (?]), the quasiparticle 
surface states on the two sides of the impenetrable wall 



are related as e B (<p) = — X ■ eJjM, where x — X l X r '■ 
Hence, if x — 1; the spontaneous currents on either side 
of the interface flow in opposite directions but with equal 
magnitude. In the case of finite transparency the sur- 
face chiral states from both sides couple into interface 
states. The absence of the total interface current on the 
scale £o is accompanied with the recovery of the symme- 
try p y — ► —p y in the spectra in Eqs.@-@ (unless $ = -k 
when there is no coupling between the surface states from 
two sides) . It is worth noting that at small transparency 
the local spectral weight of the states at x = — (or 
x = +0) is still asymmetric (see Fig. || below). 

If identical, equally oriented, massive superconductors 
on each side of the interface have opposite chiralities 
the junction can be considered as an isolating domain 
wall between two degenerate chiral phases of the bulk 
superconducting state. In this case x = ~ 1 an d spon- 
taneous surface currents in the two banks of the inter- 
face flow along the same direction. Then the recovery 
of the symmetry due to finit e tr an spa rency doesn't take 
place. Indeed, solving Eqs.(A4)-([Aq) of the Appendix 
with a model spatially constant anisotropic order param- 
eter, having opposite chiralities on the two sides of the 
junction, results in 



P,/ x 2_ 2,S nA P,/ x 2_ 

£ B \Pf)=X \\ 




\ —V cos 2 — sin (p ± v T> cos — cos ip 
z z 



(8) 



£B(P/) = -X r sgn(p y )|A*( P/ )| 



$ 



$ 



V sin 2 — cos 2ip ± vD sin — sin 2ip 



(9) 



4""(P/) 



-X r sgn(PxPy)|A{"»(p/)| 



1 - V sin 



- - — ±Vv 



■ COS <p : 



sin — sin ip 



(10) 



These branches are asymmetric with respect to the in- 
version p y — > — p y . 

The surface current, the vector potential, and the in- 
duced magnetic field at an impenetrable barrier separat- 
ing an r)i(pf)+ir)2(pf) (left side) and an 771 (p/) — ii] 2 (p/) 
(right side) superconductor are given in Figure 0. They 
decay in the depth of the superconductor on a scale given 
by the penetration depth. For the p-wave case there is 
an agreement of our results with those in Ref. |32|. If 
instead the right superconductor would have been calcu- 
lated with the same parity as the left one the currents, 
the vector potential, and the magnetic field would be re- 
versed in the right half-space. 

As was shown in the preceding subsection, the order 
parameters vary significantly from their bulk forms in 



the vicinity of the interface. This raises a question of the 
applicability of the approach used just above and based 
on a simple model of spatially constant anisotropic or- 
der parameters. We now report the results taking into 
account spatial variations of the order parameters. As 
will be seen, they display a good qualitative agreement 
with Eqs.|(p|)-(|lO|). Quantitative distinctions, however, 
take placeQ. 

The surface pair-breaking modifies the low-energy chi- 
ral bound states in a simple way. Spatial profiles of 
the order parameters near impenetrable surfaces can be 
taken into account analytically in the expressions for the 
energies, by replacing the order parameters at the surface 
by effective ones. In general, we have ^1,2(2:) 7^ and the 
expressions for the effective surface order parameters are 



quite cumbersome. We get a simple answer if the surface 
influence retains <fii{x) = 02 (a;) = 0. The last condition 
holds, in particular, for the p-wave pairing (see Fig. FH) 
and we represent here the results for this simplest case 
(the other cases are considered in the Appendix). 

Taking advantage of the approach of Refs. R,p^ (see 
also Appendix), we obtain the bound state energies in 
close vicinity of the nodes of Af(p/): 



£ s(P/) = ±X S g n ( V *MI A 21. e /7<>/' ±0 )l 



(11) 



Here and below the upper (the lower) sign corresponds 
to the right (the left) half-space; Sip is the deviation of 
ip from the direction of a node A^(p/) = (that is at 



ip = or (p — -k for the p-wave). The effective surface 
order parameters are defined as 



|Atf,e//(p/,±0)| 



±oo 



Tjttl f \ A ](Pf. x ')\ dx ' 

Ai(p f ,x)\e ° dx 



±oo 



Tj^j f \&j(Pf,x')\dx' 

e ° dx 



Eqs.(yJ), (|12|) should be compared with Eq.(|6|), 
izing it to spatially dependent order parameters 



(12) 
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FIG. 2. Spatial dependence of the surface current j y , the vector potential A y and the surface magnetic field B z . The 
penetration depth is set to 20 £o and the unit of the current is eV/NfA},. 



The transformation of the spectrum of surface bound into account if (f>i(x) = (see Appendix). For instance, 
states into one of interface states due to finite trans- for the symmetric junction (x = 1) we S e t the following 
parency of the tunnel junction, i.e. for T> <^ 1, can be result in vicinity of a momentum direction P/,0j f° r which 



described analytically taking the surface pair-breaking ertivt ) 



= 0: 



x>=o 



4(P/) = ± 



e s(P/) 



v=o 



I? cos 2 



In the case of opposite chiralities (x = — 1) we obtain 



+ OC 




^j\A\{p,,x')\dx' 

e ° dx 



1/2 



(13) 



£s(P/) = £ b(P/) 



± vDcos , 
x>=o V 2 
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FIG. 3. The local angular resolved surface density of states shown (left to right) for the p-wave, the f xy -wa,ve and the 
f x 2_ y 2-wa,ve junctions. The transparency T> = 0.01 and the phase difference is zero. The top row are all symmetric junctions 
(x = +1) while the bottom row are junctions with opposite chiralities (x — —1). The DOS is given in the interval ip — —n/2 
to <p> = 7r/2 in equally spaced steps. 



In figure || we show the numerically computed local in- 
terface DOS, using the self-consistently determined self- 
energies. The angle-resolved DOS in the left supercon- 
ductor is defined as 



N(<p, E ) = -itylmjltfoSfoO-je)]} 



(15) 



and calculated at zero phase difference with one and the 
same transparency T> = 0.01 for all momentum direc- 
tions. Here and below we use the normalization condi- 
tion g 2 = — 1 for the quasiclassical Green's function. The 
DOS in the right superconductor are for these symmetric 
junctions simply related to the DOS in the left by sub- 
stituting e r = —\ ■ e l . Qualitatively, peak positions are 



well described by Eqs.(g) and (0), obtained in the tun- 
neling limit. The fine structure of the peaks are due to 
a finite transmission. The density of states for junctions 
with bound states dispersing in pf and having e B = at 
(f = 0, i.e. for the p-wave and for the f x 2_ v 2-wave pair- 
ings, show a split DOS, by an amount \TD at e — (in 
accordance with Eqs.(||) and (||); see also ( ^3| ) and (Q). 
At finite angles ip the DOS show the predicted depen- 
dence on the relative chirality of the two superconduc- 
tors. For the x = +1 junction, since e r B v=0 = —e l B -p =0 > 
there is only a very small spectral weight in the DOS 
from the bound state s r B v induced in the left supercon- 
ductor through the junction in the tunneling limit. For 
the x = —1 junction we have e r B v=0 — e l B -p =0 . Hence, 



there is almost equal weight in the two states split around 
£b,v=o( 1 p)- The split again oc \fV (see Eq.(|i4|)). For 
small angles ((p ~ 0) the p-wave and the / x 2_ y 2-wave 
superconductors lead to qualitatively similar DOS. For 
the f X y-wave superconductors the chiral branches take on 
low energies only for momentum directions close to the 
nodes of the order parameter and the main part of the 
low-energy spectral weight belongs to glancing trajecto- 
ries. The f X y-w&ve with x = — 1 does not show splitting 
if $ = 0, in accordance with Eq.(|l0|). 



III. LOW TEMPERATURE ANOMALY OF THE 
JOSEPHSON CURRENT 

A. The Josephson current in quantum point contacts 

Consider the Josephson current across junctions with 
chiral interface states first assuming spatially constant or- 
der parameters. For a quantum point contact we can-fiad 
the current in symmetric junctions (with x — 1) asE2rO 



de% 



de% 



Here 



J = 2e£^»/(e£) = -2e £ ^f tanh^f. 

the spectrum is presumably even and the sum is taken 
over different channels. Making use of the relationship 

-4£- = — ■ — -t^- 1 which easily follows from Eqs.(||)- 

(S) for all pairings we discuss, leads to the Josephson 
current 



. e|Ar£>sin$ , e B 

1 = tanh — — 

■' 2e B 2T 



(16) 



for one channel in a quantum point contact. 

Eq. ([16]) is applicable also to contacts between isotropic 
s-wave superconductors. In this case, however, the bound 
state energies are mainly on the order of A. They can 
take low-energy values only in highly transparent junc- 
tions and for phase differences $ in a narrow vicinity of it. 
There are well-known specific features of the Josephson 
current manifested ia conventional superconductors with 
micro-constrictionsEl For the phase difference $ ~ n, 
however, sin$ in Eq.(J16|) is very small as well, preclud- 
ing the low-temperature features of the Josephson criti- 
cal current we discuss below. By contrast, under certain 
conditions chiral bound states can take low energy val- 
ues at any phase difference $. For the p-wave and the 
/ a ,2_ y 2-wave order parameters, as well as for the three- 
dimensional E lg -type of pairing, this is the case, in par- 
ticular, in tunnel junctions (T> <C 1) for quasiparticles 
with momenta aligned almost (or exactly) parallel to the 
interface normal ip w 0, tt. 

For a nonzero energy and sufficiently low temperatures, 
when T -C £b, we get from Eq.rtlq) the zero-temperature 

value of the Josephson current j = ''' "' 



"2ei 



Here 



Ao is the zero-temperature order parameter. Hence, 
the zero-temperature critical current is much greater for 



channels with eb <C |Ao|, as compared with channels, 
where Eb ~ | A, 1, for given $, V and |A |. 

If Eb <C T <C T c , which can be satisfied only by low- 
energy states, we find from Eq.(|lrj) that the Josephson 
current can vary substantially in the low-temperature 
region being inversely proportional to the temperature: 

. _ e|A | 2 I?sin$ 
3- 57 • 

To be specific, consider the p-wave order parameter 

and Eb = A^\ TZsin 2 <p + 2?cos 2 — . Here tp is the pa- 
rameter characterizing the quantum state in the chan- 
nel. The bound state energy is small, i.e. Eb <C 
for example if \ip\ <C 1 and T> <C 1. Then the zero- 
temperature current is quite large and equal to either 



A?, 



eA£ Psin$ 



1 or to 



when max 



T 

AT' 



eA? 



Dsgn 



cos -S- 



cos y 



if 



T 

AT 



< M < 



< m < 



COS-!j 



•C 1. In the last case the Josephson cur- 
rent in the tunnel junction turns out to be proportional 
to Vv, in contrast to the linear dependence on V of 
the conventional tunneling current. As we see, the low- 
temperature value of the Josephson current depends on 
the dimensionless parameter VX>|A|| cos($/2)|/T con- 
taining the transparency. This parameter appears from 
the expression for the energy of the interface bound states 
when comparing it with the temperature. At sufficiently 
low temperatures the parameter is large even if T) <C 1. 

It follows from Eqs.(||)-(p[) that the condition 
£s{Pf) <§C A(p/) is satisfied for any phase difference for 
the tunnel junctions with ip near the nodes at tpo « 0, 7r 
of the chiral bound states for the p-wave, the f x 2_ y 2-w&ve 
and the three-dimensional E\ g order parameters. For the 
f X y-wa,ve order parameter the above condition takes place 
only for glancing trajectories, which contribute negligi- 
bly to the Josephson current. For the three-dimensional 
Eiu type of pairing the nodes of the bound states, not 
coinciding with the nodes of the order parameter, lie 
at (pQ = ±7r/4, ±37r/4. For this reason the E2 U pair- 
ing (UPta) has a contribution from the low-energy chi- 
ral surface bound states to the Josephson current that 
turns out to be especially sensitive to the momentum 
dependence of the tunneling probability, T>(jjf). This 
contrasts the E\ g pairing as well as the p- and /-wave 
superconductors. The sensitivity is associated with the 
fact that the transparency of a barrier quickly dimin- 
ishes with increasing a deviation of the momentum di- 
rection from the surface normal (unless the barrier is 
sufficiently thin). Thus, the transmission can become 
quite small for momenta with ip — ±7r/4 or ±37r/4. For 
example, a model for the tunneling barrier with a uni- 
form probability distribution within an acceptance cone 
about the interface normal, and zero outside the cone, 
the effect considered is entirely absent unless the accep- 
tance cone contains (or its boundary is very close to) 
ipg = ±7r/4, ±37t/4. For a narrow Gaussian-distribution 



model the tunneling at ip^ — ±7r/4, ±37r/4 is not strietly 
zero, although exponentially small. For a sufficiently thin 
and high barrier, however, the transparency is approxi- 
mately proportional to (n-pf) 2 and then the contribution 
to the Josephson current from the low-energy chiral sur- 
face states becomes important even for the i?2 U -triplet 
pairing. In the last case the particular ^/-dependence 
of the tunneling probability is T>(pf) = VqCos 2 ip for 
two-dimensional models with a cylindrical Fermi surface 
and T>(pf) — Vq sin 2 #cos 2 ip for three-dimensional mod- 
els with a spherical Fermi surface. 



B. The critical current in classical tunnel junctions 



can noticeably weaken the low-temperature deviations of 
the Josephson current from its conventional behavior. 

The eventual result can be easily understood qualita- 
tively, if one notices that the integration of the 1/|<^|- 
tcrm (see the preceding subsection) over the interval 



I max 



T 

AT 



cos 



$ 



p c I (where (p c <§; 1) leads to 



logarithmic low-temperature dependences on the temper- 
ature or the transparency. Below we present a more care- 
ful analysis of the low-temperature features of the critical 
current. 

For a classical symmetric junction (with x = 1) be- 
tween quasi-two-dimensional superconductors we get in- 
stead of Eq. <M) : 



The low-temperature anomalous behavior of the 
Josephson current described above is quite similar to 
what was theoretically found for the case of d-wave su- 
perconductors in the presee«9-«f zero- and low-energy 
surface and interface statesc3iH3. However, the similar- 
ity takes place only in considering one channel, which 
is appropriate to quantum point contacts. In classical 
junctions, where quasiparticles have various momentum 
directions, the parameter ip appears in the expression for 
the Josephson current as an integration variable in av- 
eraging over the Fermi surface. If the bound states are 
substantially dispersive, manifesting strong dependence 
of their energy on the momentum direction, the averaging 



'^ev/iV/PMcos^sin^lA^)! 2 ^ fe B (p) 



2esM 



V 2T 



(17) 



Assuming the logarithmic functions mentioned above 
to be large, one can use a logarithmic approximation, 
considering only a part of the integral in Eq. flTfl) , which 
contains the zero-level crossing of the chiral branches 
Pffl. Then for a symmetric tunnel junction between 
quasi-two-dimensional superconductors with cylindrical 
Fermi surfaces we get from Eq.(17) 



ip +8ip a 



Y^ ev f N f V{tpo)cosp sm<P\A((po)\ 2 j 



dip- 



<P0 
P/,x>0 



»(¥>) 



tanh 



V 2T 



(18) 



¥>o 



where Eb('p) is defined by Eqs.(0) and (0), £_b(</?) is to be evaluated to linear order in Sip. Sip c is a cut-off 

parameter. The sum is taken over those po, which correspond to p/, x > 0. For the p- or the /^ ^2 -wave pairing only 
Po ~ satisfies the condition. In the case of three-dimensional superconductors with the Eiu- or the i?i ff -pairings, 

tt/2 

an additional integration J sin 2 8d8 has to be carried out in Eq.(|17[), where T>, |A| and £b, depend on both ip and 


8. For the E 2u and the E lg one should consider variations Sip from ip\ — zkir/4 and (p$ = respectively. 

The contribution to the Josephson current from low-energy bound states in the tunnel ju nction s can be also 
calculated analytically taking into account the surface pair-breaking. Indeed, from Eqs.(13), ( A14 ) we easily get 

-4£- = — ■ — ^- , where the effective order parameter |A(p/,0)| near the nodes of the bound states is defined in 

the right half space 



|A( P/ ,0)| |v/,,(p/)| 



■ —. x 2 ( P/) J \ A (Pf, x 2)\sin(ii(p f ,X2)-4>(pf : x 2 ))dx., 



dxi 



(19) 



and analogously in the left half space. Hence, spatial de- 
pendence of the order parameters introduces in this case 
the only modification in Eq.([18|): the order parameter 
|A((^)| should be replaced there by the effective surface 
order parameter |A(<p)|. 

For temperatures \fV(pf$) | A(p/,o) 1 1 cos($/2) | < 
T <C AbSipc < T c one can put tanh (£b{^p) /2T) ss 1 near 
the upper limit of the integral in Eq.rtTl) and use Eqs.(P) 

and (0) for Eg x ~ y (<p),e S g (#,</?) ex 5<p. Then we find a 
logarithmic temperature dependence of the integral near 
the lower limit (po in Eq. (p~q) . Qualitatively, the tempera- 
ture is a lower cut of the integral, but for fixing numerical 
factors we checked its asymptotic behavior numerically. 
Thus, at sufficiently low temperatures, when contribu- 
tions to the current from quasiparticles with energies 
the order of |A(p/)| already take their zero-temperature 
value, the Josephson current associated with the low- 
energy part of the chiral bound states and described by 
Eq.(|l8|) can still grow logarithmically with the decreas- 
ing temperature. As a result, the total Josephson cur- 
rent can be represented for the low-temperature interval 
v /Z>(p/,o)|A(p / , )||cos($/2)| « T « T c as 



j x (T) = aj c . ln 



bAi 
T 



sin$ 



(20) 



where j C:0 = LeN f v f V A b {LP = 2/3, L f * 2 -y 



■4 



« 0.45, 

Jxy,S,t 



LS-v = 2/5, L s = 1/9, V- = 2/21), A 6 = A^ 
a and b are constants the order of unity. In our analyti- 
cal results here and below the transparency is taken for 
high and thin potential barriers as it was represented 
in the end of the preceding subsection. Further, disre- 
garding surface pair-breaking, we calculate analytically 
a total Josephson current and find at low temperatures 
after the comparison with Eq.© a? = 3/2, bP « 0.69, 
ate-v 3 « 2.23, b f * 2 -y 2 « 0.44, a s = 3/2, b s « 0.23, 
a* = 3/4\/2, b l rs 0.95. In neglecting not only a sur- 
face pair-breaking (that is a spatial dependence of the 
order parameter), but the presence of any surface bound 
states at all, one can use spatially constant (bulk) values 
of Green's functions for the calculation of the Joseph- 
son current. This oversimplified approach leads for the 
symmetric tunnel junction to the Josephson critical cur- 
rent jefi- As the argument of the logarithmic function in 
Eq. ( P0[ ) is supposed to be quite large (for the logarithmic 
approximation to be valid) , the Josephson current j x (T) 
can noticeably exceed j c> q. 

For temperatures T -C ^/VqAi,] cos($/2)| the lower cut 
of the integral is associated with the transparency of the 
barrier and we find from Eq. (|18| ) 



j x (0) = ajcfiln 



/©o|cos($/2)| 



sin<I> 



(21) 



The logarithmic function in Eq.(pT|) is supposed to take 
large values as well. Disregarding surface pair breaking, 
we find dP sw 1.21, d f * 2 -y 2 « 0.78, d s sw 1.21, d* « 1.88. 



Based on Eqs.(|8|)-([L0|), one can easily get the Joseph- 
son current in junctions with opposite chiralities. We do 
not present here the corresponding analytical results in 
detail, since there is no unconventional low-temperature 
increase in the critical current in the case x = —1. In dis- 
regarding surface pair breaking the Josephson current in 
tunnel junction at T <C T c in the case x — ™ 1 1S written 
as 



uj cfi sin $ 



where u p — 9/8, u 



f*v = 



5/3, 



,1, 



1.95, 



(22) 
9/8, 



u ss 0.142. The "reference" current j c fi is taken here the 
same as in Eqs.(gC|) and (|2l| ) in order to make clear that 
the Josephson current in the case of opposite chiralities is 
of the conventional order of magnitude without any log- 
arithmic enhancement. The smaller numerical factor u' 
is associated with the lower transparency for ip$ = ±7r/4 

as compared with ip ' * ~" =0. Surface chiral bound 
states in the whole range of subgap energies form the 
current in the particular case, even at low temperatures. 
One can show that the low-energy chiral quasiparticle 
states lead to a current which is reduced compared to 
Eq.(p2[) by the small cut-off parameter 6ip c and it does 
not matter whether the surface pair breaking is taken 
into account or not. 

Another reason for choosing the "reference" current 
j C fi in the case x = — 1 same as for x = 1, is that in 
neglecting both the surface pair breaking and the pres- 
ence of Andreev surface bound states, and simply us- 
ing the bulk expressions for the Green's functions, the 
tunnel Josephson critical current across the domain wall 
(x = —1) vanishes. This in sharp contrast with Eq.(E2h. 
For a vanishing Josephson current one needs, for exam- 
ple, different projections of the orbital angular momen- 
tum, i.e. different values of l z , of the Cooper pairs in the 
left and the right superconductors together with a con- 
servation of l z in the tunneling process. The current van- 
ishes, for instance, if l z is along the interface normal and 
the interface is symmetric with respect to the rotations 
around the normaoa. As we demonstrate in Eq. (|22|) , if 
l z is along an axis parallel to the interface, the interface 
influence results in interface bound states and the pos- 
sibility of a nonzero tunneling of Cooper pairs. Eq.d22) 
applies to tunnel junctions. Effects of-Jaigh transparency 
can also lead to finite critical currentQ. 

On account of the bound states, the Green's function 
depends on the distance from the surface even when dis- 
regarding a spatial dependence of the order parameter. 
There is an important qualitative difference between the 
surface and the bulk values of the propagators as a func- 
tion of momentum direction. In particular, the angle <p is 
an important part of the phase not only for the chiral or- 
der parameter, but for the bulk value of the off-diagonal 
components of the Green's function as well. Since we con- 
sider the z-axis to be parallel to the surface, the surface 
breaks the conservation of the projection of the angular 



10 



momentum. For this reason the angle ip drops out from 
the phase in the surface value of the off-diagonal compo- 
nents. This can be obtained on account of the continuity 
of the quasiclassical Green's function, taken for incoming 
and outgoing trajectories (with the angles ip and n — tp 
respectively) at the impenetrable surface. Thus, on ac- 
count of the surface effects, the Josephson tunnel current 
across the domain wall is of the conventional value (see 
Eq.@). 



The logarithmic behavior of the critical current found 
above for classical symmetric junctions with x — T is ob- 
servable only at sufficiently low temperatures and at low 
junction transparencies when the logarithmic approxima- 
tion can be justified. This is demonstrated in Figure 0, 
where we show the result of numerical calculations for 
the temperature dependence of the critical current for 
different values of barrier transparency. 
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FIG. 4. The temperature dependence of the critical current divided by the maximal transparency T>o plotted for different 
values of T)q. The junctions are modelled with the angle dependent transparency T>(ip) = T)q exp( — 16sin 2 p). To the left we 
show two coupled p-wave superconductors and to the right two f x 2_ y 2-wa,ve ones. 



As seen for the x = +1 junctions, the term 2?ohi2?o 
is present in the tunneling limit in the dependence of 
the zero-temperature critical current on the transparency 
(see Eq.(f2f"l)). Further, the low-temperature range, where 
the critical current linearly depends on InT (see Eq.(p0|)), 
becomes well pronounced only for transparencies less 
than 10 -2 . The finite transparency efficiently cuts off 
the logarithmic temperature dependence of the critical 
current. The lower the transparency, the larger the low- 
temperature range. The characteristic temperature, be- 
low which the transparency influence becomes noticeable, 
is proportional to VP, again in agreement with our ana- 
lytical results. 

In Fig. |5] we show the analogous result for the f xy - 
wave superconductor. Since the low-energy contribution 
from the bound states is associated in this case only with 



glancing trajectories, there are no anomalies in the criti- 
cal current. 
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FIG. 5. The temperature dependence of the critical cur- 
rent calculated as in Fig. H but for the f xy -wave superconduc- 
tor. 



C. Effects of broadening of the bound states 

Broadening of the chiral low-energy bound states can 
substantially modify the low-temperature behavior of the 
critical current found above. Even a small broadening 
7 <C A;, have a profound influence on the low-energy 
parts of chiral branches. We take a small broadening 
into account in the pole-like term of the retarded quasi- 
classical Green's function, where we simply replace the 
factor l/[e — es(Pf)] vnth-.l/[e ~ £s(pf) + i"f], analo- 
gously to what is done inc3. This slightly modifies the 
analysis of Ref. [2l] for the contribution from the low- 
energy bound states to the Josephson current in tunnel 
junctions. We evaluate the current in the tunneling limit 
assuming the momentum independent broadening sub- 



stantially greater than the low-energy bound state and 
the temperature: T, ^/V^A b \ cos($/2)| < 7 < A b . The 
lower cut of the integral in Eq. (Hn) is associated in this 
case with 7 rather than the temperature or the bound 
state energy. Then we find for the Josephson current 



j x (0) = aj cfi ln ( I sin <T> 



(23) 



For spatially constant order parameters we get w p ~ 
0.61, w f * 2 -y 2 w 0.39, w s w 0.21, u>* « 0.84. 

There are various contributions to the broadening of 
the bound states. These are in particular associated 
with surface roughness, non-magnetic and magnetic im- 
purities, and with inelastic scattering. We assume here 
that non-magnetic impurities dominate the scattering 
and thus the broadening. Then we calculate the effect 
of small scattering rates on the behavior of the Joseph- 
son current at low temperatures using the usual t-matrix 
approximation and assuming an isotropic impurity po- 
tential, u. The impurity self energy is in this case 



E(e)=r 



(9(e)) 



l-a(l + (g(e)) 2 ) 



(24) 



where the scattering strength a — sin S s (8 S is the s-wave 
scattering phase-shift, < 6 S < 5) and the scattering 
rate, T, parameterize u and the impurity density n imp as 



(irNfi 



irNfU 2 



1 + {irNfuf 



1 + {TrNfuf 



(25) 



Due to the low-energy bound states, the Matsubara 
Green's function can take quite large values in the low- 
temperature region, if the broadening and the trans- 
parency are sufficiently small. However, the low-energy 
states form only a small part of the chiral branches, i.e. 
for most of the quasiparticle trajectories the energy of the 
chiral states is the order of A&(p/). For this reason the 
quasiclassical propagator for the chiral superconductor, 
averaged over the Fermi surface, (g(e)}, does not take 
large values but is of the order of unity or less. This 
differs greatly from the case of a 45°-oriented d-wave su- 
perconductor for which each trajectory has an Andreev 
bound state at zero energy. There (g(s)) can take large 
values along with-the pole-like term in the quasiclassical 
Green's functioned. 
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FIG. 6. The temperature dependence of the critical current for a junction modelled with the transparency 

T>((p) — T>oexp( — 16 sin <p) and T>o — 10 -5 . jc(T) is calculated at different scattering rates V/T c . The scattering strength 
is weak and put to a — 0.1. 



In Fig. B we show the influence of impurity effects on 
the temperature dependence of critical current. In or- 
der to get a wide range of well pronounced dependence 
of the critical current on the small broadening, we take 
the extreme tunnel limit and put T>q = 10~ 5 . One can 
see from the Fig. 0, that for observing the logarithmic 
low-temperature enhancement of j c , one needs supercon- 
ductors of high purity. 



nounced in quantum point contacts. In classical junc- 
tions the zero-temperature current acquires an additional 
logarithmic dependence on low transparency or on the 
broadening of the bound states. Under the conditions 
considered, the Josephson current through the domain 
wall does not vanish due to the bound state contribu- 
tion. 



IV. SUMMARY 

Chiral interface Andreev bound states have been ob- 
tained and studied above both analytically and numer- 
ically. We showed, that the low-energy chiral states re- 
sults in the low-temperature enhancement of the Joseph- 
son current between clean chiral superconductors in sym- 
metric tunnel junctions. The enhancement is more pro- 
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APPENDIX: ENERGIES OF THE INTERFACE BOUND STATES 



In the presence of a quasiparticle bound state SsiPf) the quasiclassical retarded propagator g has a pole at 
£ = £b{P}')- One can introduce the residue of the propagator g as 



g(p f ,r;e B {pf))= km [(e - e B (pf))g(pf,r;e) 



(Al) 
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which is finite and satisfies the same transport equation as g but completed with the relation 

[g(pf,r;e B (pf)j\ =0 , 



(A2) 



rather than the normalization condition. 

For calculating the bound state energies, the Eilenberger equation for g can be solved in terms of the following 
ansatzQ: 



f + {p f ,x;e B (Pf)) = g(pf,x;e B (pf))exp(-iT)(pf,x)) , 
f(pf,x;e B (p f )) = -g(p f ,x;s B (pf))exp(ir)(pf,x)) . 

The whole number of quasiclassical equations can then be reduced to the one scalar equation 



V /,z 



d x r)(pf,x)+e B (p f )- \A(pf,x)\coa(r)(p f ,x) - <f>(pf,x)) = , 



(A3) 
(A4) 



completed with the condition at the interfaceo 



'm,o{pf,i)-m,o(Pfi) 



sm 



i-i,i' 



sm 



' Vr,o(Pf,r) ~ Vr,o(Pf, r Y 



(A5) 



and the asymptotic conditions in the right and the left superconductors 

VfAPf) sin (Voo (Pf) - <t>oc(Pf)) sgna; > , 



±00 



(A6) 



E qs.(A4 )-(A6) are valid both for singlet superconductors and for triplet ones with a (0,0, d z ) order parameter. 
Eq. (|A5|) connects solutions of Eq.(A4) with momenta Pf,i, Pf, r of incoming quasiparticles from the left and the right 



f,i 



p of reflected ones. For specular reflection, the momentum parallel 



sides of the interface with the momenta p 

to the interface is conserved, i.e., p fl = pi t — p" f — p\ r . In the limit of impenetrable wall Eq.( |A5|) reduces to a 
continuity condition for f]o(pf) taken for incoming and outgoing momenta along a quasiparticle trajectory. 

Eq.(A4) can be easily solved for spatially independent order parameters. Thus, for a symmetric junction between 
the p-wave superconductors (with any chiralities) one gets 



?7/(p f; ) = $i+X;(7r-^)+arccos( £B *fJU 



-f,l J 



m(Pf,l) = ®l +Xlf~ arccosf B \J' 1, 

Vr{p fr ) = ^r + Xr^ + ^ccos{ Eb( PJ' 1 > j , r? r (p/, r ) = $ r + Xr(^ - f) - arccos ( £B ^/' l > 



(A7) 



Substitution of Eq.( A7) into Eq.( A5) results in Eq.(||) for the spectrum of the bound states in the case of identical 
chiralities and in Eq.(H) for opposite chiralities. Similar derivations with the other types of pairing lead to Eqs.(0)-(@) 
and Eqs.(|)-©. 

Calculations of surface and interface chiral bound state energies, taking into account the spatial profile of the order 
parameter, are carried out for the low-energy states with momenta close to nodes (or to low-energy minima) pf^ of 
the chiral branches e B (pf). Consider, for example, the p-w&ve order parameter. We note, that A^p/, x) = for two 
particular momentum directions: (pQ = and 7r. On the other hand, <Pq = and it correspond to the incoming and 
the outgoing momenta in a reflection event where A^(pf,x) changes its sign. Hence, the zero-energy surface bound 
states take place near an impenetrable surface for these particular quasiparticle trajectories. For the p-wave order 
parameter the solution of Eq.(AJ) for the zero-energy states is rjp ' = $ + 5sgn(x). At an impenetrable surface a 



chiral branch crosses the zero energy with a slope, which can be found for a spatially dependent orde r parameter. 
Indeed, considering a small deviation Sep = ip — (f^ from any of the two trajectories, we linearize Eq.(A4) and get the 

following first-order corrections to rjp 



(0). 



±00 



,(!) 



(x,Pf) = / (e B (p f )-x ril) sgn(p y x')\A p 2 (p f ,x')\ 



T 1 ^ T /|A»(p / ,x")|«te" 



dx 



(A8) 
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Here and below the upper and the lower signs correspond to half-spaces x > and x < respectively. 

At the impenetrable wall r](pf,x = 0) = r)p + r) p ' {x = 0,p/) is one and the same for incoming and outgoing 
momentum directions. Taking this into account we find Eq.(|Tlh for the bound state energies close to ip^ for p-wave 
pairing. 

Calculations for the -Ei g -case are very close to the p-wave pairing since the orbital parts of the order parameters 
contain identical dependence on p x and p y . Thus, we get in this case s s B (pf) = ±x s S n ( v xS'p)\A2 1 e jf (p/,±0)| with 
the same notations as in Eq.(pTj). 

Analogous consideration can be carried out for the E 2u triplet pairing. Then zero energy bound states occur at four 



particular momentum directions with tp$ = ±7r/4, ±37r/4 and this at any fixed value of 9. The solution of Eq.( A4) for 
ip = ifiQ is ry t = $ + 5- (1 + sgn(xxPyPz))- Corrections to r/l linear in small deviations Sip from ipQ, are described 
by the equation 

o f/ \ =F 1 £ T /|AS(P/,*")I*»" 

4 1 \x,p f ) = -r^J {£l(pf)sgn(y x )- X r{ %gn(x'Sip)\A\(p f ,x')\)e - dx' . (A9) 

X 

Then in a narrow region of the nodes of A\(pf) we obtain for an impenetrable wall e^(P/) = 
±XSgn(v K 5(/j)|A* 2 e tf{pf, ±0)| where 5>p is the deviation of <p from the direction of a node: A{(pf) = 0. 

If the equality (j)i(x) = 0, assumed above, is not satisfied, the result for the effective surface order parameter becomes 
more cumbersome. First, one should find the function r)(°>(x) satisfying the equation 

^-d xV ^(x) + \A 1 (p ffi ,x)\cos{7 1 ^(x)-M^-^) =0 (A10) 

and the asymptotic conditions r]( \±oo) = $ + (7r/2)sgn x. 

Then, introducing the notation d(x) = j](°>(x) — 4>i(x) — $, we obtain, for example, for the f x 2_ y 2-wave pairing the 
following result: 



/ 



- T^TT I l Al (P/' x ')l sinCi(x')da;' 

\A 2 (pf,x)\sin( 2 (x)e ° dx 
e f f- y2 (Pf) = X r( °sgn v„-3 — : — . (All) 

°° ~vh~\ j \ A i-(pf- x ')\ sin Ci-( x ') dx ' 

e o dx 



o 



Let now the transmission be finite, but sufficiently small. Assume <j>i(x) = 0, which takes place, for instance, for 
the p-wave pairing. In vicinities of the momentum directions p^ > where A?J b {p P r ) = ( that is ip = ip^ ) we get 

i] P — <& + 5-sgn(x) — sgn(xv x )e 1 ^(p P f Q )/\A^ b (p P r )\. Eq. flAgj) can be then generalized to the presence of nonzero 
transparency: 



±00 / 

'T. * 



\ ^^]\Al{p t ,x")\dx" 

-X r(i) sgn(v x x'M|Af(p/,x')|je - dx' . (A12) 

Analogously, near momentum directions P/Q' ^ 0T wn ich ^ib(Pfo) = ^ (that is <po — ±^,±^) we get rfl = 
$ + 5 (1 4- sgn(xxP y Pz)) — sgn(xVa:)(£g(pt )/l^2 b(P/ o)l- Then in the presence of a nonzero transparency we get 



/ 
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\ *t%\S \ A tot>*")\*»" 
X r{l) sea{y x x'Sip)\Ai(p f ,x')\ e - dx' 



(A13) 



For finding the bound state energies one should insert solutions from Eq.( A12| ) (Eq. (A13)) into the condition at 
the interface Eq.( A.5), taking account of the first non- vanishing corrections in transparency or in Sip. For a symmetric 
tunnel junction, when the order parameters on both sides coincide for every momentum direction and, in particular, 
have identical chiralities, this expansion leads to Eqs.(O) for the p-wave pairing as well as for the Ei g pairing. For 
the Ei u type of pairing we get in this case 



4*(P/) 



+ CO 



e*B(Pf) 




l^T / \ A i(Pf' x ')\ dx ' 
e ° dx 



1/2 



(A14) 



For the domain wall (x — —1) the spectrum for the p-wave pairing is described by Eq.(nJ). The same expression is 
valid as well for the E\ g pairing. For the E^ u type of pairing we get in the case of opposite chiralities 



<4(p/) = 4(p/) 



± VX>sin , 
v=o V 2 




~°° -^\S\^2(Vf,x')\dx' 

e ° dx 

o 



(A15) 
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